Anisotropic spin freezing in the S=l/2 zigzag ladder compound SrCu02 
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Using magnetic neutron scattering we characterize an unusual low temperature phase in orthorhom- 
bic SrCu02- The material contains zigzag spin ladders formed by pairs of S=l/2 chains (J = 180 
meV) coupled through a weak frustrated interaction (| J'| < 0.1J). For T < T c i ~ 5.0(4) K an elastic 
peak develops in a gapless magnetic excitation spectrum indicating spin freezing on a time scale 
8t > h/5E — 2 ■ 10 -10 s. While the frozen state has long range commensurate antiferromagnetic 
order along the chains (£ c > 200c) and a substantial correlation length (f a = 60(25)a) perpendicu- 
lar to the zigzag plane, the correlation length is only £b = 2.2(3)6 in the direction of the frustrated 
interaction. We argue that slow dynamics of stripe-like cooperative magnetic defects in tetragonal 
a c planes yield this anisotropic frozen state. 



PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 

Reduced dimensionality geometrical frustration 

H , and band formation [Q can suppress the critical tem- 
perature of magnets far below a cooperative energy scale 
such as the Curie- Weiss temperature. The resulting low 
temperature ordered phases have unconventional features 
that challenge conventional theories of magnetism. They 
include strongly reduced sublattice magnetization, ab- 
normal sensitivity to low levels of disorder, and mean- 
field like critical behavior. 

Quasi-one-dimensional magnetic dielectrics are excel- 
lent model systems in which to explore weakly ordered 
phases because a quantitative connection between the dy- 
namic properties of the one-dimensional units and criti- 
cal phenomena in the coupled system can be established 
^-0. Haldane spin-1 chains and other spin systems with 
a gap have a critical value for the inter-chain coupling 
zJj_ « A below which they remain disordered down to 
T = 0. Spin-1/2 chains on the other hand, are predicted 
to order even for vanishingly small non-frustrated inter- 
chain interactions. This is consistent with the low tem- 
perature long range order among weakly coupled spin-1/2 
chains found in S^CuOs and Ca2Cu03 ||. 

In this letter we describe a different low temperature 
phase in closely related SrCu02, which contains linear 
spin-1/2 chains assembled pairwise in an array of weakly 
interacting zigzag ladders |p[-fi5[. Though weak static 
sublattice magnetization does develop for T < 5.0(4)K 
« 2.8 • 10~ 3 J/fcs, three dimensional long range order is 
absent for T > 10~ 4 J/k B - Specifically, we find a low 
temperature correlation length of only two lattice spac- 
ings along the direction of frustrated intra zigzag ladder 
interactions. We argue that slow dynamics of stripe-like 
defects in quasi-two-dimensional antiferromagnetic layers 
prevent the system from developing long range order. 

Zig-zag ladders in SrCu02 are built from corner- 
sharing Cu-0 chains with an exchange constant J « 



181(17) meV [|13[ stacked pairwise in edge-sharing geom- 
etry. The frustrating interaction between chains proceeds 
through w 87.7° Cu-O-Cu bonds @|§ and is expected 
to be weak and ferromagnetic |J'| < 0.1 J. Field the- 
ory |l7|-|i"9fl and numeric simulations Jl9|-pl[ for a pair of 
chains coupled like this predict that weak antiferromag- 
netic inter-chain coupling induces incommensurate cor- 
relations and a gap A ~ Jexp(— a J/ J') while the zigzag 
chain should remain gapless for ferromagnetic J' ]18|,fL9| . 

SrCu02 is centered orthorhombic (space group 
Cmcm = D^T) with lattice parameters a = 3.556(2)1, 
b = 16.27(4)1, c = 3.904(2)1 @. We index wave vector 
transfer in the corresponding simple orthorhombic recip- 
rocal lattice. Nuclear Bragg reflections (h, k, I) are al- 
lowed for even h + k and even I when k — 0. While there 
is no direct information about the magnitude of inter- 
actions between zigzag ladders, the structural features 
that determine them are similar to those in Sr 2 Cu03. 
Application of quantum Chain Mean Field (CMF) the- 
ory t° Sr2Cu03 yields an estimate for inter-chain 
interactions of IT w k B T N / (1.28 8 J/k B T N ))_ w 
0.13 meV, while according to classical || CMF J±_ = 
3(k B T N ) 2 /{8JS 2 {S + l) 2 ) pa 8 • 10" 4 meV. 

The neutron scattering experiments were performed at 
the NIST Center for Neutron Research using cold and 
thermal neutron triple axis spectrometers. PG(002) re- 
flections were used for monochromator and analyzer, sup- 
plemented by Be and PG filters. Our sample is a cylin- 
drical rod (£ « 46 mm, D « 5 mm, and m = 3.875(5) g) 
grown by the Traveling Solvent Floating Zone (TSFZ) 
technique. Rocking curves about the b axis showed 
two roughly equal intensity peaks separated by 0.5° 
and each with a Full Width at Half Maximum (FWHM) 
« 0.25°. Experiments were performed with wave- vector 
transfer in the (h, 0, 1) and (h, fc, h) reciprocal lattice 
planes. We used incoherent scattering from a vanadium 
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FIG. 1. Wave-vector dependence of inelastic magnetic scat- 
tering from SrCuC>2. (a)-(c) show data along a* which is per- 
pendicular to the zigzag ladders while (d)-(f ) show data along 
the chains, (a)-(c) were collected on SPINS at T = 0.35K 
with E f = 5.1 meV and collimations 80' - 80' - 240'. (d)-(f) 
come from BT2 at T = 12.5(5)K with E f = 14.7 meV and 
collimations 60' — 60' — 80' — 120'. Increasing background in 
(f) comes from the direct beam. 

rod in the same geometry as our SrCuC>2 sample to nor- 
malize the magnetic scattering cross section. 

The denning characteristic of a quasi-one-dimensional 
spin system is an anisotropic dynamic correlation vol- 
ume, which can be probed by inelastic magnetic neutron 
scattering. Fig. 1 shows constant energy scans along 
two perpendicular directions in the (h, 0, 1) plane for en- 
ergy transfer 0.5 meV< fiu < 14 meV. Scans along the 
chain direction (Fig. 1 (d)-(f)) reveal a resolution limited 
peak centered at I = |. The data yield an upper limit 
Sq < 0.01c* on any intrinsic FWHM of the peak. For 
comparison, the FWHM of the des Cloizeaux-Pearson 
continuum for each antifcrromagnetic spin- 1/2 chain at 
Huj = 4 meV is predicted to be Sq = (2Tilj / n 2 J)c* — 
4.5 ■ 10~ 3 and hence it is unresolved in our measurement. 
Fig. 1 (a)-(c) show scans along the a direction which 
is normal to the plane of the zigzag ladders. There we 
find peaks for \h\ « | which span the better part of the 
Brillouin zone and this provides evidence for short range 
dynamic antiferromagnetic correlations perpendicular to 
zigzag ladders. The correlation anisotropy between the 
two directions probed is Sq a /Sq c > 10a* /c*. While our 
measurements do not yield an energy scale associated 
with dispersion along c* we found that modulation along 
a* exists only for Tilo < 2 meV. 

Susceptibility, heat capacity, and /iSR measurements 
indicate that static magnetic order develops below T ~ 
2K in SrCu0 2 g|. Figure 2 shows the corresponding 
elastic magnetic neutron scattering, which we found at 
Q = t± Q m where r is a reciprocal lattice vector, Q m = 

+ e, k, 5), e = 0.006(1), and k is an integer. However, 
rather than being concentrated in a resolution limited 
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FIG. 2. (a)-(c) wave-vector dependence of the net elas- 
tic magnetic intensity around Q « (0.5, 0, 0.5) at T=0.35K 
along b*, a* and c* directions respectively. The background 
was measured at T=8K. Horizontal bars show calculated 
FWHM of the instrument resolution function, (a) BT2 and 
BT9 measurements with Ef = 14.7 meV and collimations 
60' - 40' - 40' - 120'. The shaded area shows fc-independent 
"magnetic rod" intensity. The solid line shows Eq. (1) with 
optimized Fourier coefficients and multiplied by the squared 
magnetic form factor for Cu 2+ . Insert shows scan across the 
rod at Jfe = 0.35. (b), (c) SPINS data with Ef = 3.7 meV and 
collimations 80' — 40' — 240'. Solid lines are fits to resolution 
convoluted Lorentzians. 

Bragg peak, the signal has a characteristic distribution in 
reciprocal space, which shows that our sample does not 
have conventional three dimensional long-range magnetic 
order. The peak is sharpest in the c* direction, wider 
and incommensurate in the a* direction, and very broad 
along b*. Fig. 2 (a) and its inset also reveal a rod of 
elastic magnetic scattering extending along b*. Based 
on this information the frozen magnetic structure can be 
described as follows. Spins in each tetragonal a — c plane 
are aligned antiferromagnetically with a superposed long 
wavelength modulation or defect structure in the a direc- 
tion. Such planes are stacked with a correlation length 
of only a few times the lattice spacing along b. 

By fitting data in Fig. 2 to Lorentzians duly convo- 
luted with the instrumental resolution function, the mag- 
netic correlation lengths in the a — c plane were deter- 
mined to be £ a = 60(25)a and £ c > 200c . The peak 
position refined to Q m = (0.506(1), 0.00(1), 0.500(1)). 
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Note that the incommensurability along a* was repro- 
duced in several independent experiments. Moreover, we 
found an equivalent magnetic satellite from r = (202) at 
Q = (1.494(3), 0.00(1), 1.500(1)) « r - Q m . 

The extremely short correlation length along b sug- 
gests analysis of scans along this direction in terms of a 
Fourier series: 



S aa (Q) 



(i 



-L^C^cosQ-(d 3 -d,0), (1) 



where index consecutive a — c spin planes and 

Cj,r = ^/N\\)Y. m {S d]+ J{S A]l+ J/(S)\ For simplic- 
ity we have assumed that the frozen spin configuration 
is isotropic in spin space and an exponential decay of 
correlations between planes separated by distances of 2b 
and beyond. The line through the data in Fig. 2 shows 
the result of this fit. The correlation length extracted is 
£b = 2.2(3)6. The correlation parameters Cjf are listed 
in table 1. 

SrCu02 is built from copper bilayers. Equivalent spins 
within a bilayer are displaced by di do = (OSh) to 
form zigzag ladders. Bilayers are stacked in registry - 
ABABAB- where neighboring A and B type bilayers are 
separated by d2 — do = (||0). Table 1 shows that de- 
spite the proximity of planes that make up a bilayer, 
correlations between them are weak. This is experimen- 
tal evidence for an effective decoupling between zigzag 
ladder rungs. Looking beyond a bilayer we see antifer- 
romagnetic correlations between spins in shifted bilayers 
A and B (C i, C03, C 5, C07 < 0) and ferromagnetic corre- 
lations between spins in bilayers separated by multiples 
of b (C 3,Co±4,Co5,C 7 > 0)- Though they are separated 
by half the distance, nearest neighbor A-B correlations 
are significantly weaker than A-A correlations and this is 
evidence for frustrated inter-bilayer interactions. From 
the fit we also obtain a reliable value for the frozen mo- 
ment at T = 0.3 K which turns out to be only a minute 
fraction of the full moment: g < S >= 0.033(7) per Cu. 

The temperature dependence of the squared frozen mo- 
ment and the in-plane correlation parameters are shown 
in Fig. 3. The frozen staggered magnetization appears 
below T c \ = 5.0(4)K and initially increases in proportion 
to (T c i - T) 2/3 where = 0.46(12). Below an inflection 
point at T C 2 = 1.5(3)K, (S) 2 increases faster until sat- 
urating below T ~ 0.5K. Careful inspection of specific 
heat data jl^] for SrCu02 reveal anomalies close to both 
T c i and T C 2- While the peak position along the chain, 
Iq, is temperature independent, both the incommensura- 
bility and the Half Width at Half Maximum (HWHM), 
k/j = a/(jr£ a ) along a*, decrease by approximately a fac- 
tor 2 below T C 2- It is useful to compare our results to 
those for Sr2Cu03 M. The critical temperature (T/v = 5 
K) and intra-chain exchange constant (J = 220 meV) are 
similar for the two materials. Nonetheless, the moment 
in SrCuC>2 is almost twice smaller than in S^CuOa. 
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FIG. 3. Temperature dependencies of (a) the frozen sub- 
lattice magnetization squared, (b) the in plane wave vector 
(Q = (ho-O.lo)), and (c) the inverse correlation length along 
a. Data derived from h and I scans through the (i + e 0^) 
elastic magnetic peak on SPINS using Ef — 4.47meV and 
collimations 80' — 80' — 240'. The frozen moment in (a) is 
within error bars of our global average. 

There is also a qualitative distinction between ordering in 
these materials; the frozen phase in SrCu02 is incommen- 
surate along a, while Sr 2 Cu03 has commensurate order 
||. Magnetic disorder along b is also unique to SrCuO^. 
While this feature may be connected to quenched disor- 
der, possibly from oxygen non-stoichiometry in our sam- 
ple, the magnetic correlation length along b is easily an 
order of magnitude less than the defect spacing. SrCu02 
therefore falls in the category of magnets which are cither 
intrinsically disordered or display an extreme sensitivity 
to disorder. We shall argue that highly frustrated intra- 
and inter-zigzag ladder interactions are central to un- 
derstanding these differences, and the anisotropic frozen 
state in general. 

The main contribution to inter-zigzag ladder interac- 
tions along a comes from weak superexchange along Cu- 
O-O-Cu paths with two » 90° bonds. For each spin 
there are 12 such equivalent Cu-O-O-Cu paths to neigh- 
bors in the same a-c plane; four to each nearest neigh- 
bor along a, denote the corresponding exchange constant 
J' a > 0, and one to each next nearest neighbor in diag- 
onal (1,0, ±1) directions, J" a > 0. As a consequence 
we expect that J" a ~ </' a /4 corresponding to a highly 
frustrated point for this lattice. Similar interactions exist 
in the (1, 5, directions between nearest neighbors in 
adjacent a — c planes. In fact, it is evident that these 
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interactions as well as the inter-bilayer interactions all 
share the same frustrating zigzag geometry and all cancel 
at the mean-field level. In addition there may be other 
weak non-frustrated interaction favoring long range or- 
der. Experimental evidence for competing interactions in 
SrCuC>2 lies in the incommensurability of the magnetic 
peaks along the a* direction. We note that the shift of 
the magnetic peak along a* equals its HWHM not only 
at the lowest T but also as a function of T < T C 2 . This 
indicates that we are dealing not with a periodic modula- 
tion superimposed on otherwise antiferromagnetic a — c 
planes but a highly disordered or even random sequence 
of frozen stripe defects induced by these competing in- 
teractions. 

Recent theoretical studies |||2| have established that 
coupled S=l/2 chain systems can remain disordered at 
T=0 if interchain interactions are sufficiently frustrated. 
In the quasi-two-dimensional case which may be a good 
first approximation for SrCu02, disorder at T=0 could 
result from instanton topological defects created through 
quantum tunneling 22 1. For half-integer or odd integer 
spin systems the disordered ground state is predicted to 
be degenerate and gapless while even integer spin systems 
should be gapfull. Our experimental results support a 
gapless phase for S=l/2. 

At T > topological disorder is enhanced through 
thermal creation of defects, leading to a reduced correla- 
tion length. In the coupled S=l/2 chain system we ex- 
pect defects to be strongly anisotropic, taking the form of 
stripes extending along c between bands of phase shifted 
antiferromagnetic domains. For order to develop in a 
stack of a — c planes, stripe defects in neighboring planes 
must move into registry. However, as mentioned above, 
the residual interactions which favor such ordering are 
either exceedingly weak or cancel at the mean field level. 
Thus it seems plausible that pinning and/or intrinsically 
slow dynamics of such stripe defects might lead to spin 
freezing instead of long range order. 

In summary we have found an anisotropic spatially 
disordered frozen spin configuration among interacting 
zigzag spin- 1/2 ladders in SrCuC>2. The magnetic corre- 
lation length along the short direction of the zigzag ladder 
is far less than the impurity spacing. Highly frustrated 
interactions both within the zigzag spin-1/2 ladders and 
between ladders, as well as slow dynamics and pinning 
of order-destroying stripe defects are likely reasons that 
SrCu02 behaves differently from the closely related lin- 
ear chain system S^CuOa. Previous theoretical work 
on coupled S=l/2 chains has indicated the possibility of 
an intrinsic disordered phase P pl]]^ . Given our data 
it would be interesting to further explore this possibility 
with competing interactions of the specific type found in 
SrCu0 2 . 
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TABLE I. Correlations between spins in a — c planes dis- 
placed by dji — do with respect to each other. Coordinates 
are given in units of orthorhombic cell parameters. S = 0.122 
[jy] is the spacing in units of b between linear spin chains that 
form a zigzag chain. Error bars define intervals wherein x 2 
is statistically indistinguishable from its minimum value. "?" 
indicates correlations whose contributions to Eq. (1) cancel 
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1 (OS I) -0.09(6) 
-1 (|<5-§§) -0.09(5) 
±2 (i±|0) ? 

3 + -0.19(5) 
-3 (05 -If) 0.06(5) 
±4 (0±10) 0.23(3) 


5 (01 + c5±) 0.09(4) 
-5 (§<5-f§) -0.18(4) 
±6 (|±|0) ? 

7 -0.11(4) 
-7 (OS -2±) 0.06(4) 
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